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Abstract 

It is shown that a light beam in free space is representable by an integral over a vectorial angular 
spectrum that is expressed in terms of an extension matrix, which describes the vectorial nature of 
the beam. A symmetry axis of the extension matrix is identified. When it is neither perpendicular 
nor parallel to the propagation axis, we arrive at such beams that show us for the first time the 
observable evidence of the Beauregard effect. The advanced representation theory may yield any 
kinds of light beam, and the uncovered Beauregard effect would play its unique roles in applications. 
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Introduction The vectorial feature, or the po 
beam, which concerns its angular momentum 



: polarization property, of a coherent light 
[ll; l2l, [3( ; intensity distribution , and 



diffraction characteristics M, l8[], has become more and more important in diverse areas of 

— ^-rwa — '—"Oft— 1 

trapping and manipulation [12], [13[, and dark-field imaging [14|. For a plane wave or a nearly 



plane wave, we have the Jones vector [15( to describe the vectorial feature of the electric 
or magnetic field. But for a bound beam, we have not yet had such a clear concept for its 
vectorial feature. 

In this Letter, we advance a representation theory of free-space light beams in terms of 
an extension matrix acting on a two-form angular spectrum. The extension matrix in this 
theory plays the role of describing the vectorial feature of a bound beam. A symmetry axis of 
the extension matrix is identified due to the transversality property of the electromagnetic 
wave in free space. What is surprising is that when the symmetry axis of the extension 
matrix is neither perpendicular nor parallel to the propagation axis of the beam, we arrive 
at such beams that show us for the first time [id] the observable evidence of the Beauregard 
effect. By Beauregard effect we mean a novel phenomenon, surmised more than 40 years 
ago by Beauregard [17| , that a circularly polarized beam of zero transverse wave vector can 
be deflected in a transverse direction. It is shown that the transverse deflection due to the 
Beauregard effect can be as large as the order of the beam waist under specific conditions. 
The uncovered Beauregard effect would open a new area of applications. 

General theory We consider three-dimensional light beams that propagate in free space. 
Let us first write out, in the Cartesian coordinate system, the integral expression for the 
vectorial electric field of a beam that may propagate in positive x direction, 

E(x) = / / ^f^A(k y , k z ) exp(zk ■ x), (1) 

J Jkl+k 2 z <k 2 Zl1 

where the time dependence exp(— iut) is implied and omitted. As we know, any element 
of the angular spectrum represented by its wave vector k = ( k x k y k z ) T has only two 
independent polarization states, though its vectorial amplitude A = ( A x A y A Z ) T has 
three components, where the superscript T means transpose. Denoting respectively by s and 
p the two orthogonal linear polarization states, we have A = A s + A p = A s s + A p p, where 
A s and A p are the complex amplitudes of the s and p polarization states, respectively, s and 
p are their respective unit vectors. Letting s = s x e x + s y e y + s z e z and p = p x e x +p y e y +p z e z , 



where e X} e y , and e z are the unit vectors in the directions of the Cartesian coordinates, Sj 



and pj (j = x, y, z) are real numbers, we have 



A = PA, 



(2) 



where 



A 



A s 

An 



(3) 



161 ] of the angular spectrum which is, or can be regarded as, the 



is the two-form amplitude 

Jones vector [l^] associated with the angular spectrum, and 
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(4) 



extends the two-form amplitude A onto the three-component vectorial amplitude A and is 
thus referred to as the extension matrix. It is clear that it is the extension matrix rather 
than the two-form amplitude that describes the vectorial nature of the beam. The mutual 
orthogonality between the unit vectors s and p and the wave vector k leads to 

2 I 2 I 2 i 

2 i 2 i 2 i 
Px+Py+Pz = !, 

SxPx + SyPy + S z p z = 0, (5) 

^x^x ~t~ ^ySy ~t~ ^z^z 0} 

kxPx + kyPy ~l~ kzPz = 0. 

We have represented a free-space light beam by Eqs. ([I])-©. But there are only the above 
five equations to determine the six real numbers Sj and pj in the extension matrix. That is 
to say, there is one freedom to be chosen. This can be done by specifying the symmetry axis 
of the extension matrix. Two independent types of representation will be discussed below. 

Without loss of generality, we consider the following two-form amplitude, 



A 



h 



A = IA, 



(6) 
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when' 1 = [ | describes the polarization state of the angular spectrum and is assumed to 

satisfy the normalization condition v l = 1, the superscript f means transpose conjugate, and 
A is the amplitude distribution of the angular spectrum (ADAS). For the sake of simplicity, 
we will consider the following Gaussian distribution function [id ]. 



A=(^) 1/2 exp 



W 2 W 2 2 



(7) 



where w y = wq/ cos 6*o, w z = wq, u>o is half the width of the beam at waist, and k y o = k sin 8q. 
This ADAS indicates that the axis of propagation, represented by its principal wave vector 
k = ( k x o k y0 ) T , is perpendicular to the z axis, where k x0 = kcos9 . For the same 
reason, only beams satisfying paraxial approximation 18j A6 = <C 1 will be considered, 



where A8 is half the divergence angle of the beam. Under this condition, Eq. ([I]) can be 
rewritten as [3] 

E(x) = — / / A(ky, k z ) exp(ik • x)dk y dk z . (8) 

J— oo J -co 

Representation of p z = and uniformly polarized beams When p z = 0, we have for the 
other elements of the extension matrix that solve Eq. (J5j), 

k x k z 
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" ~ k(k 2 x + k 2 y y/ 2 ' 

(k 2 + k 2 ) 1 / 2 
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X 



(k 2 x + k 2 y y/ 2 ' 
Py = W+WP' 

Denoting k r = k r e r = k x e x + k y e y in a cylindrical coordinate system, where k x = k r cos 9, 
k y = k r sin 6, and e r is the unit vector in the radial direction, it is clear that p = p x e x +p y e y = 
eg is the unit vector in the azimuthal direction. Furthermore, letting s r = s x e x + s y e y , it 
is found that s r = — -ye r is in the radial direction and that s z — ^. All these facts show 
that the symmetry axis of the extension matrix (Tj0) is the z axis and thus is perpendicular 
to the axis of propagation. Because the ADAS in Eq. ((7j) is appreciable only in a small 
region in which \k y — k y o\ ~ \k z \ < kA6 <C k, k z /k in the extension matrix can be regarded 



as a small number in comparison with unity in view of the integral (jHJ). To the zeroth-order 
approximation, we have for the extension matrix, 



' -sin O ^ 
cos 6> 



V 



o 



(9) 



/ 



The electric field corresponding to extension matrix (Q is 

E T (x) = [l x e z + I 2 (e y cos9 - e^sinflo)] 



— J J Aexp(zk ■ x.)dk y dk z 



(10) 



which is transverse with respect to the propagation axis ko and is uniformly polarized. The 
polarisation state is represented by that of the angular spectrum, ( Higher-order correc- 
tions will be included if higher-order approximations are considered [4|. When appropriate 
ADAS's are selected, the known kinds of uniformly polarized beam will be obtained, includ- 



ing 



19 



;he fundamental Gaussian beam, Hermite- Gaussian beams, Laguerre-Gaussian beams 



201 ]. and the Bessel-Gaussian beams 
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Representation of s x = 0, non-uniformly polarized beams, and the Beauregard effect We 
have in this representation the following solution to Eq. (jSJ), 

A', 
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Denoting k r = k r e r = k y e y + k z e z in another cylindrical coordinate system, where k y 



k r cos(p, and k z = k r sinip, we find that s 



Sy^y Sz^z 



e,„ is the unit vector in the 



azimuthal direction. In addition, letting p r = p y e y + p z e z , it is apparent that p r 



is in the radial direction and that p x = — %■. That is to say, the extension matrix in this 
representation is symmetric with respect to the x axis. Therefore, 9$ here is the angle 
between the symmetry axis of the extension matrix and the propagation axis. It will be 



shown in the following that different values of 6q correspond to different beams with different 
polarization properties and intensity distributions. All the beams in this representation are 
spatially nonuniformly polarized. 



The case of 8 = and axially symmetric polarization The symmetry axis of the 
extension matrix is parallel to the axis of propagation in this case. We use the cylindrical 
coordinate system by defining x = xe x + r, where r = ye y + ze z , y = r cos 0, and z = r sin 0. 
The ADAS in Eq. (jTJ) is expressed in the cylindrical coordinate system as 



Substituting Eqs. 
expansion, 



), (HU), and (H2J) into Eq. 



exp (ip cos ip) = i m J m (p) exp(imip) 

m=—oo 

we obtain for the electric field 



(12) 

and with the help the following 

(13) 



where 



Ei(x) = z(Zie^ + l 2 e r )rxi{r, x) exp(ikx) + l 2 e x Ef x (x), 



r n 4-ur 



(14) 
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w(x) = w ( 1 + i 
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J m 's are the Bessel functions of the first kind, Jo and I\ are the modified Bessel functions 
of the first kind. In deriving Eq. f[T4"|) . we have made (i) the paraxial approximation [3] 
k x w k — if in the exponential factor exp(ik • x), (ii) and the zeroth-order approximation 



p r = %e r ~ e r in the extension matrix. 

Eq. ( [HP describes beams of axially symmetric polarization. The first term on the right 
side is the transverse component and is of the zeroth order. The second term is the longitudi- 
nal component. It is of the first order, ~ k r /k, and is thus much smaller than the transverse 



component 18|, |23(. Neglecting the small longitudinal component, the beam is dark on the 
principal axis r = 0. Locally, it is elliptically polarized the same as the angular spectrum. 
Other kinds of axially symmetrically polarized beams will be obtained when appropriate 
ADAS's are selected, includin g th e azimuthally and radially polarized Bessel- Gaussian and 



Laguerre-Gaussian beams 



24 



25 



26|. 



The case of \9 \ < A6> and the Beauregard effect If \9 \ is not equal to zero but is 
very small satisfying l^ol < AO, we have approximately w y ~ wo, so that the ADAS in Eq. 
([7]) becomes 

A = A exp(k y0 wlk y ), (15) 



where 



A W ° 

A = —= exp 

'71 



W 
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Ko ) ^P ( -=^k 



u2 



With the help of expansion (Tl3|) and the Graf formula 27] . we obtain for the beam's electric 
field after substituting Eqs. ®, JTTJ, and (JTSD into Eq. ©, 

E 2 (x) = E 2T (x) + he x E^), (16) 

where 



E 2T (x) 



A 2 
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[i(he^, + l 2 e r )r + (l 1 e z + l 2 e y )k y0 wl 

xx 2 ('~f,x)exp(ikx), 
f°° K 

— exp(ikx) / -^A 2 Jo('yk r )k r dk ri 
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A exp y-^K i ■ 
(r 2 — ky WQ — 2ik y0 wlr cos^) 1 / 2 . 



In deriving Eq. ( fl6l) . we have also made the aforementioned two approximations. 

The first term on the right side of Eq. (fT6|) is the transverse component, and the second 
term is the longitudinal component, which is also of the first order as before. It is interesting 



to note that the transverse component consists of two parts having different polarization 
symmetries. One is of axially symmetric polarization and is dark on the propagation axis 
r = 0. The other is uniformly polarized. But locally both of them have the same polarization 
state as that of the angular spectrum, I. It is the interference between those two parts that 
yields the Beauregard effect as will be shown below. 

Neglecting the small longitudinal component, the intensity distribution of the beam is 
found to be 

I 2 « |E 2T | 2 = (y 2 + z 2 + 2ak y0 w 2 z + k 2 y0 w*)\ X2 \ 2 , (17) 

where a = i(lll2 — hi 2) is the polarization ellipticity of the angular spectrum. The first factor 
on the right side of Eq. (|17[) shows that the intensity is minimum at a point r = —e z ak y0 WQ 
on the z axis. The displacement of this point from the propagation axis is 

z = —ak y0 WQ = —akwQsmOo. (18) 

For angular spectra of circular polarizations a — ±1, the minimum intensity is equal to 
zero. The single-sided depression of the beam intensity on the z axis renders the beam 
centroid deflected transversely in the z direction unless a or #0 is equal to zero. According 
to Eq. (fT8l) . the deflection will be flipped to the opposite side of the propagation axis by 
changing either the sign of a or the sign of 9q. Since the expectation value of the transverse 
wave vector k z of the beam vanishes, the transverse deflection of the beam centroid in the z 
direction is nothing but the Beauregard effect [16;]. To the best of my knowledge, this is the 
first time to demonstrate the observable evidence of the Beauregard effect. As an example, 
we show in Fig. [I] the normalized intensity distribution of the beam on plane x — 0, where 
a — 1, and #0 = A# = 10~ 3 rad. The intensity in Fig. [I] is numerically calculated directly 
from Eq. ( 1T61) and is well approximated by Eq. (1171) as numerical calculations show. The 
null point is located at zq ~ — wq = — 159A on the z axis. As a result, the beam centroid is 
deflected in the positive z direction by (z) ~ 102A as can be obtained by using formula (24) 



of Ref. 



16), which is as large as the order of wq. 



The case of \8 \ > A6 The electric field distribution of the beam in this case cannot 
be approximately expressed so clearly as have been done above. But the propagation char- 
acteristics of a paraxial beam can be studied numerically by use of the integral (jSJ). For a 
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500 



FIG. 1: Normlized intensity distribution on plane x = for a = 1 in representation s x = 0, where 
9q = A8 = rad, the y and z coordinates are in units of the wavelength A. 

given #0; the beam deflection in the transverse direction due to the Beauregard effect can 
be calculated by using formula (24) of Ref. Q. 

In conclusion, the representation theory advanced here may produce any kinds of free- 
space light beam that solve the Maxwell equations. The Beauregard effect uncovered here 



would be closely related to the Imbert-Fedorov effect [16j and open a new area of applications, 
"t also has vital implications to the counterpart of the quantum-mechanical matter waves 
281 ] . In order to explore the Beauregard effect in experiments, one should be able to control 
separately the propagation axis of the beam, the symmetry axis of the extension matrix, as 
well as the polarization ellipticity of the angular spectrum. 

This work was supported in part by the National Natural Science Foundation of China 
(60377025), Science and Technology Commission of Shanghai Municipal (04JC14036), and 
the Shanghai Leading Academic Discipline Program (TO 104). 
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